Photonic heterostructures with Levy-type disorder: statistics of coherent transmission 
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We study the electromagnetic transmission T through one-dimensional (ID) photonic heterostruc- 
tures whose random layer thicknesses follow a long-tailed distribution -Levy-type distribution. 
Based on recent predictions made for ID coherent transport with Levy-type disorder, we show 
numerically that for a system of length L (i) the average (— InT) oc L" for < a < 1, while 
(— InT) oc L for 1 < a < 2, a being the exponent of the power-law decay of the layer-thickness 
probability distribution; and (ii) the transmission distribution P{T) is independent of the angle of 
incidence and frequency of the electromagnetic wave, but it is fully determined by the values of a 
and (InT). Additionally we have found and numerically verified that (T) oc L~" with < a < 1. 

PACS numbers: 03.65.Nk, 42.25.Dd, 72.15.Rn 
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Random processes characterized by density probabil- 
ities with a long tail (Levy-type processes) have been 
found and studied in very different phenomena and fields 
such as biology, economy, and physics. One of the main 
features of a Levy-type density distribution p{l) is the 
slovif decay of its tail. More precisely, for large I, 



P{1) 



1 



ll+o 



(1) 



vifith < a < 2. Note that the second moment diverges 
for all a and if < a < 1 also the first moment diverge. 
This kind of distributions are also knovifn as a-stable dis- 
tributions [1]. A vifindow on new optical materials vifhich 
allow for the experimental study of Levy flights in an 
outstanding controllable way was recently opened with 
the construction of the so-called Levy glass [2]: titanium 
dioxide particles are suspended in a matrix made of glass 
microspheres. The distribution of the microsphere di- 
ameters is properly chosen in order that light can travel 
performing Levy flights within the microspheres. The 
diameter distribution is characterized by the exponent a 
of the power-law decay of its tail; it was found [2] that 
when < a < 1 the transport is supperdiffusive, while for 
a = 2 the normal diffusive transport is recovered. This 
experimental investigation has motivated several theo- 
retical works on the effects of the presence of Levy-type 
processes on different transport quantities in one dimen- 
sion, as well as in higher dimensional systems [3-6, 8]. 

On the other hand, coherent electron transport 
through one-dimensional (ID) quantum wires with Levy- 
type disorder was studied in Ref. [6]. It was found that 
for the dimensionless conductance, or transmission, T: 

(i) the average (over different disorder realizations) of 
the logarithm of the transmission behaves as 



(-InT) 



for 
for 



0< a < 1 
1 < a < 2 



(2) 



and 



(ii) the distribution of transmission P{T) is fully deter- 
mined by the exponent a and the ensemble average 

(InT). 

We point out that although for 1 < a < 2, the average 
(InT) depends linearly on T, as in the standard Ander- 
son localization problem, it is interesting to remark that 
the statistical properties of T are not those predicted 
by the standard scaling approach to localization, in par- 
ticular by the Dorokhov-Mello-Pereyra-Kumar (DMPK) 
equation [7]. That is, for 1 < a < 2 the transmission fluc- 
tuations are larger than those considered in the DMPK 
equation. Thus, the standard statistical properties of T 
are recovered for a > 2. 

Having as a reference various analogies between elec- 
tron, light, and matter-wave transport [9-13] one may 
expect the statements (i) and (ii) to be also valid for 
ID optical systems, even though in the latter case addi- 
tional parameters such as incidence angle and frequency 
come into play. A systematic investigation of the statis- 
tical properties of coherent transmission for a ID ana- 
log of Levy glasses is not available in the literature and 
the above expectations have not been verified until now. 
Therefore, in this paper we undertake this task by study- 
ing the transmission T through ID photonic heterostruc- 
tures with Levy-type spacing disorder. Moreover, engi- 
neering the disorder in a photonic heterostructure might 
be a less complex task than in an electronic system; so 
we hope that this work stimulate future photonic exper- 
iments. 

The heterostructures that we shall study consist of an 
alternating sequence of layers of materials A and B hav- 
ing refractive indexes ua and ns, respectively. See Fig. 1. 
The corresponding thicknesses Ia and Ib are chosen as 
random numbers drawn from a Levy-type distribution 
characterized by the exponent a of its power-law tail. In 
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FIG. 1: (Color online) Sketch of a heterostructure of length L 
with random layer thicknesses Ij and refractive indexes rij^^gy 
The distribution of the thicknesses follows a Levy-type distri- 
bution. 



this work wc only consider Levy-type distributions sup- 
ported along the positive semiaxis and focus on the case 
< a < 1, where large fluctuations of the transmission 
produce the most interesting effects. The total length L 
of the heterostructure is then given by i = ^ ^j . With- 
out lost of generality we consider that the heterostructure 
grows in the z-direction. 

ID heterostructures with Levy-type thickness disorder, 
as defined above, can be produced with porous silicon. 
Layers of silicon with different porosities (i.e. different 
refractive indexes) and varying thicknesses can be ob- 
tained by HF electrochemical etching by modulating the 
value of the current density during the anodization pro- 
cess [14]. Thus we fix ua = 1.4 and ub = 2.4 in our 
calculations since they correspond to experimentally ac- 
cessible porous silicon refractive index values [14, 15]. 

We compute the transmission T through our ID het- 
erostructure using the transfer matrix formalism de- 
scribed in Ref. [16]. We consider an electromagnetic wave 
with frequency w that strikes the first layer of the het- 
erostructure (embedded in air) making an angle 9 with 
respect to the z-axis, then propagates inside the struc- 
ture composed of N layers, and finally escapes through 
the layer on the opposite side. Without loss of generality, 
in the following we specialize on TE modes. The transfer 
matrix of the scattering process can be written as [16] 



M = Mo.ArMnetMi n 
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where M^ 
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and kjz is the component of the wave vector along the 
z— direction in the jth layer given by kjz = kjCos{6j) 
with kj ~ {to/cjrij [uj equals ua or ns)- Oj is related to 
dj-i through Snell's law: sin{9j)/ sin{9j-i) = nj-i/rij. 
Above, Mi^o [Mq^at] is the transfer matrix for the inter- 
face between the first [A^th] layer and air, while Mh„t is 
the transfer matrix of the heterostructure. Finally, 



In the following we shall study the statistical properties 
of T, in particular its full distribution. Concerning the 
numerical simulations along this work, the statistics is 
collected over an ensemble of different disorder realiza- 
tions. Notice that since the thicknesses li are drawn from 
a Levy- type distribution, Eq. (1), for a fixed length L the 
number of layers composing a heterostructure might vary 
strongly from sample to sample. Along the paper, uj is 
given in units of the reference frequency ojq ~ 2'n:c/XQ 
where Ag could be chosen to provide suitable experimen- 
tal conditions. 

As we have mentioned, the theoretical results in this 
paper are based on the analysis presented in [6]. Here 
we only reproduce the main result of that work: the dis- 
tribution for transmission -P^(r) with ^ = (InT) . The 
average quantities that we study, such as (T) and (In T) 
[Eq. (2)] are derived from the following expression for the 
distribution 



PdT) 



Ps{a,i,x)iT)qaA{x)dx 



(5) 



for < a < 1, where qa^c is the probability den- 
sity function of the Levy-type distribution supported 
in the positive semiaxis, s{a,£_,x) = ^/{2x"Ia), la = 
1/2 Jp°° x^'^Qa^idx, and 



Ps{T) 



s ^/^ exp(-s/4) 
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(6) 
where j/o = arcosh(2/T — 1). 

Wc start by analyzing the average transmission (T). 
From Eq. (5), we find that 



(T) oc L- 



(7) 



T = |M2 



(4) 



for < a < 1. We have verified numerically this result. 
In Fig. 2 we present the average transmission for different 
values of 9 and tu ioi a — 1/4 and a = 1/2. As we can 
see, the agreement is very good in all cases. This result, 
Eq. (7), could be contrasted with the exponential decay 
of (r) with L, for standard ID disordered systems, and 
with (T) ex l/L, for quasi ID systems in the normal 
diffusive transport regime. 

We now study the average (hiT). In electronic trans- 
port this quantity is of relevance since it gives information 
about the localization length of the disordered system. In 
Fig. 3 we present different plots of the average (hiT) as a 
function of L for a < 1 [Fig. 3(a-b)] and a > 1 [Fig. 3(c- 
d)]. For a < 1 we observe a clear behavior of the form 
(—InT) oc i"; while for a > 1 we see that (—InT) is 
simply proportional to L, as for standard ID disordered 
systems. Therefore, all curves in Fig. 3 are well described 
by Eq. (2). In addition, notice that while the curves of 
(InT) vs L strongly depend on the incidence angle 9 for 
all a, see Fig. 3(a) and 3(c), the dependence on a; is lost 
for a > 1, see Fig. 3(b) and 3(d). 

Next, we analyze the full distribution of the transmis- 
sion. It is clear from Fig. 3 that by choosing the appro- 
priate combination of 0, uj and L one can fix the value 
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FIG. 2: Average transmission (T) as a func- 
tion of L (symbols) for ID heterostructures with 
Levy- type spacing disorder characterized by a. 
In (a) [(b)] three values of 9 [uj] were considered 
for ljJ — 0.25 [S = 0]. Dashed curves are fittings 
of the numerical data accordingly to Eq. (7). 
Each symbol was calculated using 10* ensemble 
realizations. 
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FIG. 3: The average — (InT) as a function of 
L (symbols) for ID heterostructures with Levy- 
type spacing disorder characterized by a. In 
(a,c) [(b,d)] three values of 9 [u] were considered 
for ui = 0.25 [9 = Q\. The dashed curves are fit- 
tings of the data with Eq. (2). Each symbol was 
calculated using 10 ensemble realizations. 



panel RED BLUE GREEN BLACK 



(a) 


(0,0.1) (7r/5,0.5) 


(^/4,10) (3^/7,1) 


(b) 


(0,10) (7r/5,l) 


(7r/6,0.1) (7r/3,0.5) 


(c) 


(0,10) (7r/5,5) 


(27r/5,l) (7r/6,0.1) 


(d) 


(0,2) (7r/4,10) 


(^/5,1) (^/3,5) 



TABLE I: Values of {9,u]/u)q) used for the histograms (labeled 
by their color) in Fig. 4. 



of the average (InT) for a given a. Then, in Fig. 4(a- 
b) we show probability distribution functions P{T) for 
a ^ 1/2 and 3/4 for (-InT) = 2. Since for smaller 
values of (InT), P{T) is concentrated close to T = 0, 
in Fig. 4(c-d) we present P(lnT) for (-InT) w 8. No- 
tice that each panel in this figure contains 4 histograms 
with different combinations of 6 and uj. However, all 



histograms fall one on top of the other; that is, P{T) 
and P(lnT) are completely determined by a and (InT). 
Moreover, the black dashed lines are the corresponding 
theoretical predictions for P{T) and P(lnr) for the spe- 
cific combinations of a and (InT) we used. Evidently, 
the correspondence between theory and numerics is ex- 
cellent affirming the equivalence between quantum and 
electromagnetic ID disordered systems. 

Finally, we want to stress that even though we have 
used arbitrary units for the length L of the heterostruc- 
tures, our results may be experimentally confirmed by 
properly choosing the frequency ui = l-KajX. For in- 
stance, since w/ojo — Aq/A, if Aq = 2L a typical experi- 
ment in the visible range with A ~ 500nm and L ~ 10/x 
[14, 15] will set the ratio w/wq to 40. Then, Q can be 
tuned to specify a desired value of (T) or (In T) . 

In summary, for ID heterostructures whose layer thick- 
nesses follow a Levy-type distribution, characterized by 



a =1/2 
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FIG. 4: (Color online) The probability distri- 
bution functions P{T) (a-b) and P(lnr) (c-d) 
for a = 1/2 and 3/4 (histograms). Each panel 
contains 4 histograms with values of 9 and to 
given in Table I. The histograms in (a-b) [(c-d)] 
are characterized by - (InT) « 2 [- (InT) ?a 8]. 
Black dashed curves are the corresponding the- 
oretical predictions obtained from Eq. (5). Each 
histogram was calculated using 3 x 10* ensemble 
realizations. 



a power-law decay, we have shown that: (—InT) oc L" 
for < a < 1; and once a and (InT) are fixed the dis- 
tribution of transmission P{T) is invariant with respect 
to the system parameters {0,uj). Also, we have found 
that (T) ex i~", in contrast to the exponential decay in 
standard disordered ID systems. We have verified that 
our results are unaffected by considering TM modes and 
different refraction index contrasts ua/tlb- 

We hope that our results serve to deepen our under- 
standing of transport properties when Levy- type disorder 
is present and motivate further experimental investiga- 



tion on its unconventional effects, as those shown in this 
work. 
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